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ABSTRACT 

WMAP's detection of high electron-scattering optical depth r e suggests substantial 
star formation at high redshift z ~ 17±5. On the other hand, the recovered as ~ 0.84± 
0.04 disfavors a cluster Sunyaev-Zeldovich (SZ) origin for the observed small-scale- 
CMB fluctuation excess, which generally requires as ~ 1.1. Here we consider the effects 
of high-redshift star formation on the CMB. We derive a fairly model-independent 
relation between r e and the number of ionizing photons emitted per baryon iV 7 , and 
use this to calibrate the amount of high-redshift supernova activity. The resulting 
supernova remnants Compton cool against the CMB creating a Compton-y distortion 
y ~ few x 10~ 6 within observational bounds. However they also create small-scale SZ 
fluctuations, which could be comparable with SZ fluctuations from unresolved galaxy 
clusters. This raises the exciting possibility that we have already detected signatures 
of the first stars not just once, but twice, in the CMB. 



1 INTRODUCTION 



The recent detection of high electron-scattering optical 
depth r = 0.17 ± 0.04 by the Wilkinson Microwave 
Anisotropy Probe (WMAP) suggests a reionization redshift 
z r = 17 ±5 (Kogut et al 2003; Spergel et al 2003), provid- 
ing good evidence for significant star formation (SF) at high 
redshift z. WMAP combined with other large-scale structure 
data also supports a ACDM cosmology with power-spectrum 
normalization as = 0.84 ± 0.04. 

This power-spectrum normalization is discrepant from 
that inferred from the CMB-fluctuation excess at small 
scales (Mason et al 2002; Dawson et al 2002), if this excess 
is attributed to the Sunyaev-Zeldovich (SZ) effect from un- 
resolved groups and clusters (Bond et al 2002; Komatsu & 
Seljak 2002; Goldstein et al 2002). These observations re- 
quire a 8 (n b h/0.035) - 29 = 1.04 ±0.12 at the 95% confidence 
level (Komatsu & Seljak 2002). 

It has been argued that galactic winds could give rise 
to a detectable SZ effect (Majumdar, Nath & Chiba 2001). 
Here we argue that the stellar activity required to photoion- 
ize the Universe at z r ~ 20 injects a considerable amount of 
energy into the IGM, which is then transferred to the CMB 
due to the efficiency of Compton cooling at these high red- 
shifts. Although the resulting mean Compton-?/ distortion is 
consistent with the experimental upper limit, there may be 
detectable angular fluctuations in the y distortion. We show, 
in fact, that for reasonable reionization parameters the fluc- 
tuation amplitude from high-z SF may be comparable to 
that from galaxy clusters. If so, then the above-mentioned 
discrepancy in the power-spectrum normalization may be 
resolved. 

In the next Section we argue that supernova remnants 



at z > 10 cool by Compton heating of the CMB and discuss 
the energetics of this process. In Section 3, we derive a rela- 
tion between the measured optical depth r e and the number 
of ionizing photons required to reionize the Universe. We 
then show that this number of ionizing photons is propor- 
tional to the energy injected into the IGM by supernovae, 
and thus the energy transferred to the CMB. In Section 4 
we discuss angular fluctuations in the y distortion and show 
that they may be comparable at small scales to those from 
unresolved clusters. 

In all numerical estimates, we assume a ACDM cos- 
mology given by the best fits to the WMAP data: 
(fi m ,fi A , fi 6 , h, a 8 ) = (0.27, 0.73, 0.044, 0.7, 0.84). 



2 HOW DOES THE SUPERNOVA REMNANT 
COOL? 

At redshifts z > 7, galactic winds powered by multiple (> 
10 ) supernovae (SN) or an energetic quasar jet are cooled 
primarily by Compton cooling from the CMB (Tegmark, 
Silk & Evrard 1993; Voit 1996; Madau, Ferrara & Rees 
2001). Less powerful winds result in cooler remnants where 
radiative losses could potentially be important. However, at 
z ~ 10 — 20, the wind from even a single SN will lose a 
substantial fraction of its energy to the CMB, as we show 
below. 

Zero- met allicity stars should be supermassive, M* > 
IOOM0, due to the thermodynamics of molecular-hydrogen 
(H2) cooling (Abel, Bryan & Norman 2000; Bromm, Coppi 
& Larson 2002). Furthermore, pair- instability SN from such 
Very Massive Stars (VMSs) should have explosion energies 
~ 100 times more powerful than conventional type II SN, 
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E V ms ~ 10 53 erg (Heger & Woosley 2002). An extreme but 
plausible version of zero-metallicity SF in low-mass halos 
T V i r < 10 4 K prevalent at high z is 'one star per halo', where 
internal UV photodissociation of H2 by the first star in that 
halo halts all further gas cooling and SF (Omukai & Nishi 
1999; Glover & Brandt 2001). Simulations show the VMS 
quickly photoevaporates all the gas within the shallow halo 
potential well within a sound-crossing time (M. Norman, pri- 
vate communication). In addition, the VMS photoionizes a 
region around the halo up to R ~ 7O(M V ms/1OOM ) 1/3 kpc 
comoving, assuming that each baryon in the VMS can ion- 
ize ~ 10 s HI atoms (Bromm, Kudrizki & Loeb 2001). Thus, 
the SN remnant (SNR) expands into a pre-ionized region at 
roughly the mean IGM gas density. During the adiabatic 
Sedov phase, R = 'y {Et 2 / pigm) 1 ^ 5 , where 70 = 1.17. The 
remnant is no longer adiabatic and begins to Compton cool 
when t ~ tc, where the Compton cooling time is, 

tc = Zm^AoTaT^y 1 = 1.4 x 10 7 [(1 + z)/20]- 4 yr, (1) 

independent of temperature and density. The (proper) size 
of the remnant at this point, when it quickly loses most of 
its energy, is 

R = 2.2 (iJvA/s/10 53 erg) 1/5 [(1 + z)/20]- 11/5 kpc (2) 

in physical units. The angular scale is 9 = R/cIa = 
0.9" (E V Ms/10 53 erg) 1/5 [(l+z)/20] -11/5 (which corresponds 
to I = tv/0 = 7.6 x 10 5 ), beyond the reach of present- 
day CMB interferometers. Thus, SNRs are effectively point 
sources, unless many SN explode together in the same 
galaxy, and/or SN bubbles from clustered halos overlap (see 
below) . 

Most of the mass and energy of the remnant is in the 
dense postshock shell, which is at p s h e ii ~ 4piGM- At t = tc, 
we can compute the temperature behind the shock front 
from the Sedov- Taylor solution, v s = 0.47 o (i?/piGMi 3 ) 1/ ' 5 , 
and assuming a strong shock T c = 3v 2 fj,m p /16kB- We thus 
obtain the ratio of Compton and isobaric radiative cooling 
time t r ad = 2.5fc s T/(nA(T)) at t = tc as 



tc 



= 0.4 



Evms 
10 53 erg 



)0.4 „ 
V 20 



a; 



(3) 



where A(T) = A23 10 ergs cm , and A23 ~ 1 for low- 
met allicity gas with T ~ 10 5 — 10 7 K. Thus, roughly a third 
of the SNR energy is lost to Compton cooling. 
The electron-ion equilibration timescale 



id = 10 yr 



15 



/ T \ 3/ 2 

U0 6 K, 



(4) 



(where 5 is the overdensity of the postshock shell) is signif- 
icantly shorter than the Compton cooling time at all red- 
shifts, so there is no problem in quickly transferring the 
shock energy from protons to electrons. 

A perhaps more likely scenario is one where many 
stars M Mot ~ 10 7 (/«/0.1)(/ 6 /0.1)(AW/10 9 ) M (where 
fb ~ flb/Q m is the baryon fraction, and /» is the frac- 
tion of baryons which fragment to form stars) form to- 
gether in rarer, more massive halos T v i r > 10 4 K where 
atomic cooling allows much higher gas densities and more 



efficient SF (Oh & Haiman 2002). The massive-star evolu- 
tion timescale is i» ~ 3 x 10 6 yr <C tc- Thus, if SF takes 
place in a starburst mode, the explosions are essentially si- 
multaneous, and Etot ~ NsnEsn- Then, an extremely ener- 
getic wind powers a much hotter bubble, and from equation 
©, trad/tc oc 2> 1 and radiative cooling is entirely 

negligible. For instance, if /* ~ 10% of the baryons in a 
Mdm ~ 10 9 M Q halo fragment to form VMSs, tc ~ trad/40). 
In principle, radiative losses could be significant in the dense 
ISM of these larger halos (since photoevaporation does not 
take place in these deeper potential wells); however in prac- 
tice most simulations (e.g. Mac-Low & Ferrara (1999)) find 
that for such low-mass systems, the SN bubbles quickly 
'blow out' (particularly in disks) and vent most of their en- 
ergy and hot gas into the surrounding IGM. Hereafter we 
shall encapsulate this uncertainty as e ~ 0.3 — 1, the av- 
erage fraction of the explosion energy lost to the CMB via 
Compton cooling. If stars form in clusters in higher-mass 
halos rather than singly in low-mass halos we expect this 
efficiency to be high, e > 0.8. 

The spatial distribution of SF does not affect 
our estimate of the mean Compton-y distortion (y = 
ksThot / (m e c 2 )Thot)- more clustered SF results in higher T e 
but lower r^ot- However, it does of course affect the strength 
of SZ fluctuations. We now turn to these issues. 



3 THERMAL SUNYAEV-ZELDOVICH 
EFFECTS 

3.1 SZ flux from Individual Supernovae 

The SZ flux from an individual SNR is 



2fr 3 T 2 

S„ = -7§-^g{x) I dQ\ATu(e)\ 



2fc|r 3 



h 2 c 2 
= 1.8 x 10 



k B T e N, 

g{ x ) — 
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where g(x) = a; 4 e 3: [a;coth(a;/2) — 4]/(e x - l) 2 is the spectral 
function, x = hv/kT-,, T 7 = 2.7K is the CMB temperature, 
and N e is the total number of hot electrons at temperature 
T e . In the second line we have used kbT e N e » Evms- The 
energy of the remnant is a function of time, EvMs(t) ~ 
£ ; VMS,oexp(— t/tc) (in the regime where Compton cooling 
off the CMB dominates). The flux from an individual SNR 
is well beyond threshold for any realistic experiment; only 
a very large number of SN (> 10 s ) going off simultaneously 
within a star cluster will be detectable at the ~ mjy level. 
Thus, SN bubbles cannot be identified and removed from 
SZ maps; unresolved SN will create both a mean Compton- 
y distortion and temperature fluctuations, which we now 
calculate. 



3.2 Mean Compton y distortion 

We first use the observed optical depth r e to derive a lower 
limit to the number of ionizing photons iV 7 emitted per 
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baryon. The dominant contribution to t oc (1 + z) comes 
from high z where the recombination time t rec oc (1 + z)~ 3 
is short, and recombinations are the rate-limiting step to- 
ward achieving reionization. The filling factor of HII regions 
is Qhii ~ t rec /t ion ps N-[/(aBn e (z)Cii(z)), where JV 7 is 
the rate at which ionizing photons are emitted per baryon 
(in units of s _1 ), Cn = (n 2 e ) / (n e ) 2 is the clumping factor 
of ionized regions (e.g., Madau, Haardt & Rees (1999)). 
The clumping factor increases with time as structure forma- 
tion proceeds; it declines sharply at high z and is C ps 2 at 
z = 20, compared to C « 30 at z = 10 (Haiman, Abel & 
Madau 2001). More sophisticated considerations (Miralda- 
Escude, Haehnelt & Rees 2000) take into account the den- 
sity dependence of reionization, but apply primarily near 
the epoch of overlap, Qn — * 1, when overdense regions are 
ionized. This has little impact on our estimates. Most of the 
mass and the optical depth comes from regions close to the 
mean density. 

The electron-scattering optical depth is given by: 

f , dt f N-, \ ca T N-, 

r e = car / dz—n e (z)mm 1, = -^—.(6) 

J dz \ a B n e (z)Cu J a B Cn 

Due to the cancellation of the electron density, this expres- 
sion is independent of the redshift of reionization, and the 
evolution of the comoving emissivity N-y(z) with redshift; 
it allows us to directly relate r e and iV 7 . The only redshift 
dependence lies in the effective clumping factor C'n, which 
increases if reionization takes place at late times. The second 
equality breaks down if overlap Qu — > 1 is achieved at high z 
and N 1 /{aBn e {z)Cn(z)) > 1 (i.e., recombinations no longer 
balance ionizations); in using the expression we would then 
underestimate N 7 , which would only imply an even larger 
emissivity. The high optical depth r e = 0.17 ± 0.08 (2a) 
(Kogut et al 2003; Spergel et al 2003) detected by WMAP 
therefore implies that 

N^ GM = 17 ±8 (T/10 4 ) ~ ' 7 (C/j/4) (7) 

ionizing photons were emitted per baryon, where T is the 
mass- weighted temperature of the reionized IGM (the y~ () - 7 
factor arises from the temperature dependence of the re- 
combination coefficient). Consistency with WMAP requires 
more SF if reionization took place at lower redshift, due to 
the increase in gas clumping at late times. 

Since only a fraction / osc of ionizing photons escape 
from their host halo due to photoelectric absorption, the 
actual total number of ionizing photons produced is larger, 
iV*°* = N^ GM f^l . In addition, we only care about those 
photons emitted at z > 6, when tc < tu (where tu is the 
Hubble time) and Compton cooling is most efficient. Since 
T e (z < 6) ?a 0.05, we have r e (z > 6) ps 0.12; therefore 
0.12/0.17 ~ 0.7 of the photons are emitted at z > 6. Thus: 

N* ot (z > 6) « 25 (/csc/0.3)- 1 (r e (z > 6)/0.12) (C/j/4) . (8) 

Estimates for the escape fraction span / csc ~ 10~ 2 — 1, but 
if the earliest stars reside in low-mass halos with T vir < 10 4 
K, the gas in such halos is quickly photoionized and driven 
out in a photoevaporating wind. If so, / csc ~ few x 0.1 to 

fese ~ I- 

How much SF and energy production is associated with 



N* ct ? We consider first VMSs, supported as the source of 
reionization perhaps by elemental-abundance evidence from 
low-metallicity halo stars (Oh et al 2001) and theoretical 
modelling (Cen 2002; Wyithe & Loeb 2002). Bromm, Ku- 
drizki & Loeb (2001) find that for 300M Q < M» < 1000M., 
the luminosity per solar mass is approximately constant; 
for M* ~ IOOMq, it falls by a factor of 2. Our estimates 
are thus independent of IMF details. For 1 ionizing pho- 
ton per baryon in the universe, /» ~ 10~ 5 baryons have 
to be processed into VMSs; thus, N^ ot = 25 corresponds 
to /» ~ 2.5 x 10~ 4 . A ~ IOOMq pair-instability SN releases 
10 53 erg (Heger & Woosley 2002), or E b ~ 0.5MeV 
per baryon processed into the VMS. The total energy release 
per baryon is therefore: 

E c = ef^Eb = 100 (e/0.8) (N* ot /25) eV, (9) 

where e is the fraction of the thermal energy which is lost to 
the CMB. A possible caveat is if a large fraction of the mass 
in the first stars went into VMSs with M, > 260Mq , which 
may collapse directly to black holes without exploding as 
SN (Heger & Woosley 2002). 

The fraction of baryons processed into VMSs /* ~ 2.5 x 
10~ 4 (iV 7 /25) implies an IGM metallicity Z~6x 10~ 3 Z Q , 
assuming uniform enrichment (since ~half the VMS mass 
is thought to end up as metals). This is consistent with the 
observed metallicity of the Lye* forest at z = 3 of Z ps 
10~ 2 5 Zq, which is not observed to evolve strongly at higher 
z (Songaila 2001). Thus, the metals seen in the Lya forest 
may well have been injected at very high z by Pop III stars. 
No trace of the entropy injection associated with the metal- 
polluting winds would remain, due to the high efficiency of 
Compton cooling. 

Our derived ionizing-photon:energy:metal ratios would 
also hold for normal stellar populations (rather than VMSs), 
which produce roughly the same amount of SN energy and 
metals per ionizing photon. The arguments are also roughly 
independent of IMF, as the massive stars that emit ionizing 
photons also eventually explode as SN. 

We now compute the Compton-y parameter associated 
with this energy injection. For simplicity, we assume that all 
of the energy is injected at some redshift Zi. The actual red- 
shift evolution introduces at most a factor ~ 2 uncertainty 
(see expresion below). The y parameter is then given by 

y = [ ° dtn e (t)E c , e- (t -^ /tc (10) 

E c 

ps n e (zi)a T ctc(zi) %r 

m e c 2 

= 3.6 x 10~ 6 ((1 + Zi/ 15 )" 1 (£ c /100eV) , 

where we have moved the electron density outside the in- 
tegrand, n e (t) ~const, since the density does not change 
significantly on the timescale over which the gas Comp- 
ton cools. In the RJ limit, (AT /T) = -2y = 7 x 10~ 6 . 
The y distortion is less than the COBE FIRAS constraint, 
y < 1.5 x 10~ 5 (Fixsen et al 1996), as it should be. Such a 
y distortion could in principle be detected by future instru- 
ments (Fixsen & Mather 2002). In addition, a low- frequency 
distortion due to free-free emission from ionized halos should 
also be detectable (Oh 1999). 
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We pause here for a simple order-of-magnitude check. 
Let the total amount of energy per baryon injected through 
Compton cooling into the CMB be E c . If this takes place at 
some median redshift Zi, this introduces an energy density 
perturbation of the CMB AU 1 ~ n b E c ~ 6.8 x 10~ 2 ([1 + 
z]/15) 3 (£ c /100eV)eVcm" 3 . The CMB energy density is 
C/ 7 = 1.3 x 10 4 [(l + z)/15] 4 eVcnT~ 



resulting in a tempera- 



ture perturbation, 



AT 



1 AU 



5.2 x 10 



V 15 J UOOeW 1 ' 



roughly consistent with our previous estimate, from 
(AT/T) = — 2y. Why is the mean y distortion due to non- 
gravitational heating by high-z SN competitive with that 
from galaxy clusters today? By integrating over the Press- 



Schechter mass function and assuming T ga3 



T 



find that the mean mass-weighted gas temperature today 
is (T) = 0.7 keV. However the Compton cooling time in 
clusters is tc ~ 150£h, so only e ~ 1/150 of that energy is 
extracted. Since y oc E c (l + z)~ 4 , we find that the y dis- 
tortion due to clusters is ~ (0.7keV/150)/0.1keV X 15 ~ 1 
times the distortion due to high-z SN. 



4 SZ FLUCTUATIONS 

Angular SZ fluctuations can be induced by Poisson fluctu- 
ations in the number density of sources, as well as by clus- 
tering of the underlying mass distribution. Poisson fluctua- 
tions are the dominant source of SZ fluctuations for galaxy 
clusters, but they are negligibly small for high-z SN for the 
following reason: like high-z halos, clusters are ~ 2 — 3<r 
fluctuations at the epoch at which they form and contain 
roughly the same fraction of collapsed mass; however, they 
are more massive by ~ 6 orders of magnitude and hence 
have a much lower space density (in addition, the comov- 
ing volume in the local universe is smaller). Thus, for the 
same Compton-?/ parameter, the Poisson contribution to an- 
gular CMB fluctuations will be negligible compared to that 
from SZ clusters. We have verified this by direct numerical 
calculation. 

We thus turn to the CMB fluctuations induced by 
clustering of high-z SN. We suppose for simplicity that 
stars form only in halos where atomic cooling can oper- 
ate, T V i r > 10 4 K, where some constant fraction /» of the 
baryonic mass fragments to form stars. Both of these are 
normalized to produce the same total fraction of baryons 
processed in VMSs f*, g iobal = 2.5 x 10~ 4 (iV 7 /25). We use 
Press-Schechter theory to calculate the abundance of halos. 
A hot bubble around each source has a total flux S oc Esn 
as given by equation @, and lasts for a Compton cooling 
time tc- The size of the hot bubble is given by equation 
©. The finite bubble size damps the power spectrum on 
scales below the bubble size. For simplicity we shall assume 
yi = j/ exp[— (l/lc) 2 ], where yi is the Fourier transform of 
the y profile of the bubble, l c = n/6 c and 8 C is the angular 
size of the bubble when most of the Compton cooling takes 
place. If the SF efficiency is independent of halo mass then 
?/o = RMhaio, where the normalization constant R oc ef* is 
determined from the condition that: 



y = I dz(dV/dzdQ) 



dM{dn/dM)y {M,: 



(12) 



subject of course to the condition that f co u > fvMS and 
ef, < 1. 

In reionized regions, gas accretion is suppressed in halos 
with T vir < T min x 2.5 x 10 5 K (or u cir ~ 50km s _1 , Thoul & 
Weinberg (1996)); lower-mass halos are thus unlikely to be 
able to form stars. This boosts the clustering bias of SF sys- 
tems as reionization proceeds, which increases the strength 
of SZ anisotropies. To keep our analysis general, we con- 
servatively only require T v t r > 10 4 K, but then show how 
increasing the Jeans mass would boost the clustering bias 
thereby enhancing CMB fluctuations. 

The Compton y power spectrum due to clustering of 
sources is given by: 



Ci(y) = / dz 



dV 
dzdQ. 



P(k 



d M (z) ' 



(13) 



dn 



dM — b(M,z) yi (M,z) 



where P(k) is the linear power spectrum, du = g(a(1 + z) 
is the comoving angular diameter distance, and b(M, z) is 
the linear bias factor (Mo & White 1996). We have used 
the Limber approximation k = Ijdu which is valid for small 
angles. Note that Ci{AT/T) = 4:Ci{y)in the RJ limit. The 
results are shown in Figure Q for two cases: (A) a standard 
"best-estimate" case with y — 3.6 x 10 -6 and clustering bias 
associated with T V i T > 2 x 10 4 K halos; and (B) a maximal 
case with y = 10~ (consistent with the current uncertainty 
in r e , C11, and / eS c), the largest value allowed by the COBE 
constraint y < 1.5 x 10 -5 ), and clustering bias associated 
with T V i r > 10 5 K halos. Also shown are the cluster-induced 
power spectra for = 0.84 ± 0.08 (2a), computed as in 
Cooray (2000). Although the "best-estimate" reionization 
signal lies below the cluster signal, with current uncertainties 
they could plausibly be comparable. The shape of the power 
spectra are fairly well constrained, but their amplitude is 
uncertain by ~ 1 — 2 orders of magnitude, as we discuss 
below. 

Roughly speaking, the CMB power spectrum is Ci ~ 
y 2 wi, where y is the mean Compton y parameter from equa- 
tion llll and wi oc l n (if P(k) oc k n ) is the flux-weighted halo 
angular power spectrum. The flatness of wil' 2 at high / is 
because P(k) oc k~ 2 at these wavenumber. For halos with 
T v ir > 10 4 K, the rapid increase in bias tend to cancel the 
decrease in the growth factor at high z, and the halo corre- 
lation function and power spectrum b{M 1 )b(M 2 )D(z) 2 P(k) 
(D(z) is the linear-theory growth factor) do not evolve 
strongly with redshift. We see this in Fig. [5] where we plot 
[b(M(T c ),z)D(z)] 2 , and b{M{T c ),z) is the mass-weighted 
bias, 



b(M c ,z) = 



dM—Mb(M,z 
dM y 




dM d JL M , (14) 



which corresponds to the flux-weighted bias since we assume 
S oc M. This is likely a minimal estimate of the bias since 
the SF efficiency (and hence the thermal SZ flux) is likely 
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Figure 1. The small-scale power spectrum of the CMB. Our 
standard case A for thermal SZ from reionization assumes y = 
3.6 X 10 -6 and clustering bias associated with T v i r = 2 X 10 4 K 
halos, while a maximal case B assumes y = 10 -5 , and clustering 
bias associated with T v i r = 2 X 10 4 K halos. A minimal case (not 
shown) which assumes a lower y or lower clustering bias would be 
down at the 10~ 14 — 10 — 13 level. Also shown are the power spectra 
due to physics at the surface of last scatter, and the thermal-SZ 
effect for clusters for its = 0.84±0.08 (2<r) (dashed lines are curves 
for erg = 0.76,0.92). 




z 



Figure 2. The evolution of the normalization factor for the halo 
power spectrum, [D(z)b(M, z)] 2 , where D(z) is the growth factor 
and b(M, z) is the mass- weighted bias from cquation ll4l The halo 
power spectrum for halos where T v i r > 10 4 K (which can cool 
by atomic cooling) hardly evolves with redshift; the increase in 
bias cancels the decrease in the growth factor. In addition, as 
reionization proceeds, halos with T v j r < 2.5 X 10 5 K are unable 
to accrete gas; hence the clustering bias of SF halos will increase 
further. The thermal-SZ power spectra in Fig 1 scale directly with 
this factor. 



to increase with the depth of the potential well. As reion- 
ization proceeds, the actual bias interpolates between the 
two curves, since accretion is suppressed in halos forming 
in reionized regions with T„j r < 2.5 x 10 5 K; it approaches 
the upper curve as Qn — » 1. Since we are probing scales 
on order of or smaller than the halo correlation length, 
r ~ few Mpc comoving, it is reasonable to expect pro- 
jected halo density (and hence flux) enhancements of order 
1(1 + 1)wi/(2ty) ~ few. 

Overall, our primary uncertainties in the predicted am- 
plitude are due to uncertainties in the mean y parameter 
that arise from the uncertainties in r e , C'n, and / eS c dis- 
cussed above. There is then an additional uncertainty of 
~few introduced by the range of halo bias factors illustrated 
in Fig.H 



5 CONCLUSIONS 

We have pointed out that WMAP's large electron-scattering 
optical depth r e implies that SZ fluctuations from high-z SF 
could be considerable. As an interesting secondary result, we 
derive a relation between r e and A r 7 , the number of ioniz- 
ing photons emitted per baryon. We use this to calibrate 
the amount of SN activity, and thereby obtain the expected 
Compton-y distortion, y ~ few X 10 -6 . Fluctuations in the 
Compton-j/ parameter could be detectable and may well 



account for the the small-scale CMB-fiuctuation excess at 
small angular scales. If so, small-scale CMB measurements 
are not a reliable independent measure of ag. If the small- 
scale CMB anisotropics are due to clusters alone, they will be 
resolved by forthcoming high-sensitivity and high-resolution 
SZ surveys. On the other hand, if high-z SF contributes sig- 
nificantly, there will be a substantial unresolved component, 
since the extremely faint flux from individual halos is unde- 
tectable. A large amount of high redshift SN activity also 
produces X-rays (Oh 2001), with interesting consequences 
for reionization. 

If a high-z origin of the observed small-scale CMB fluc- 
tuations is confirmed, CMB maps may then be used to study 
the topology of reionization, perhaps by cross-correlating 
with future 21cm tomographic maps of neutral hydrogen 
at high z (Tozzi et al 2000). Here we have focused exclu- 
sively on thermal-SZ fluctuations, which induce a Compton- 
y distortion to the CMB frequency spectrum and can thus 
be distinguished from "genuine" temperature fluctuations 
with multifrequency CMB measurements. However, high-z 
SF may also induce temperature fluctuations by scattering 
from reionized regions with coherent large-scale peculiar ve- 
locities, as we detail in an accompanying paper (Cooray et 
al., in preparation). 

Given the uncertainties in high-z SF discussed above, 
we can make predictions for small-scale y fluctuations with 
roughly an order-of-magnitude level of uncertainty in the 
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CMB-fluctuation amplitude, and thus cannot at this point 
conclusively attribute observed small-scale CMB-fluctuation 
excesses to high-z star formation. Nonetheless, this interpre- 
tation of the excess is certainly plausible. If it is correct, 
then the CMB experimentalists have achieved a remark- 
able triumph: not only have they fulfilled a decade-old quest 
to measure cosmological parameters with exquisite and un- 
precedented precision, they have detected signatures of the 
very first generation of star formation not just once, but 
twice. 
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